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Abstract 

We discuss two heuristic ideas concerning the spectrum of a Laplacian, and we 
give theorems and conjectures from the realms of manifolds, graphs and fractals 
that validate these heuristics. The first heuristic concerns Laplacians that do not 
have discrete spectra: here we discuss the notion of "spectral mass" closely related 
to the "integrated density of states", an average of the diagonal of the kernel of 
the spectral projection operator, and show that this can serve as a substitute for 
the eigenvalue counting function. The second heuristic is an "asymptotic splitting 
law" that describes the proportions of the spectrum that transforms according to 
the irreducible representations of a finite group that acts as a symmetry group 
of the Laplacian. For this to be valid we require the existence of a fundamental 
domain with relatively small boundary. We also give a version in the case that the 
symmetry groups is a compact Lie group. Many of our results are reformulations 
of known results, and some are merely conjectures, but there is something to be 
gained by looking at them together with a unified perspective. 
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1 Introduction 



What is "spectral asymptotics" if not the study of the asymptotic behavior as A — t- oo 
of the eigenvalue counting function 



N{X) = #{A, < A} (1.1) 

of some operator? For this to make sense the operator must have a discrete spectrum. 
In this paper we restrict attention to operators that may be called "Laplacians" , but the 
classical theory [RS] is not restricted to Laplacians, and many of the ideas in this paper 
can be extended at least to Schrodinger operators (Laplacian plus potential) and perhaps 
beyond. But what if the operator does not have a discrete spectrum? Must we give up 
on the idea of spectral asymptotics? Not at all! 

We will discuss a notion of "spectral mass" M(A) that can replace A^(A). To motivate 
the definition recall the famous Weyl asymptotic law 



N{X) = cn vol(fi)A"/2 + o(A"/2) as A ^ oo (1.2) 

when the operator is —A on a compact domain Q in M" (or more generally a Riemannian 
manifold of dimenson n) with smooth boundary with Dirichlet or Newmann boundary 
conditions. (Here c„ is an exphcit dimensional constant.) We divide fll.2p by vol(f2), and 
note that N{X)/ vol(f2) is equal to the average value of Kx{x, x) on fl, where Kx is the 
kernel of the spectral projection operator Ex onto the [0, A] portion of the spectrum, 

Kx{x,y) = ^j{x)^j{x) (1.3) 

A,<A 

where {fj} is an orthonormal basis of eigenfunctions with associated eigenvalues {Aj}. 
So in this case we define the spectral mass function 



^(^) = — TT?^ / Kxix,x)dfiix), (1.4) 
vol(i2) Jq 

{n is Lebesgue measure or the Riemannian mesure on Q). Our simple observation above 
is that 

M(A) = (1.5) 
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and so Weyl's asymptotic law says 



M(A) = c„A"/2 + o(A"/2). 



(1.6) 



Now consider a Laplacian so that —A is a nonnegative self-adjoint operator on 
L^(f2,/i) where the measure is not finite. The spectral theorem still supplies the spectral 
projection operators Ex, and in many cases these are given by integration against a kernel 
Kx{x, y), but of course (11.31) is not valid. Our first heuristic idea is to try to define M(A) 
as some sort of average value of K\{x, x) and to study its asymptotic behavior. Note that 
we do not expect (11.61) to hold, and indeed in many cases the asymptotic behavior will 
be as A —i- or some other natural value. Also, in some cases we may have to settle for 
upper and lower (M"'"(A), M~(A)) spectral masses if we can't prove that a true average 
exists. This is clearly a heuristic idea because it is necessary to decide on the meaning 
of "average" in each context. 

We should point out that the same idea has been used in quantum mechanics under 
the name integrated density of states. The usual context is the Schrodinger operator 
—A + V{x) on Euclidean space for a potential V{x) that is periodic [BSh], almost 
periodic [Sh],[PS], or random ([KM] and references therein). The most common definition 
is to take 



where fi^ is a sequence of regular domains increasing to M", and Nk is the eigenvalue 
counting function for —A + V{x) on Qk with suitable boundary conditions (in the peri- 
odic case it is typical to take Qk to be cubes with periodic boundary conditions). The 
advantage of this definition is that it allows for numerical approximations, it avoids deal- 
ing with the spectral projections for operators that do not have discrete spectrum, and 
it implies monotonicity with respect to the operator. It is shown in [Sh] that for almost 
periodic potentials the definition agrees with our definition (2.5), where the average is 
the Bohr mean of almost periodic functions. The main point of this paper is to export 
this idea to many new contexts, and this is the justification for introducing the term 
"spectral mass" in place of "integrated density of states." 

In Section 2 we will elaborate on the idea of a spectral mass, establishing its ele- 
mentary properties (positivity and monotonicity), discuss its relationship to the heat 



N{X) = lim 



1 



iVfc(A), 
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kernel, and present a conjecture related to the minmax characterization of eigenvalues 
that would allow simple estimates of spectral masses. 

In Section 3 we discuss examples where the Laplacian is homogeneous, so that K\{x, x) 
is constant and there is no problem defining the average. 

For Euclidean space, we obtain (11. 6p exactly (no error term), while for the lattice 
we obtain (11.61) as A — )• (the spectrum is bounded). For hyperbolic space we obtain 
(1.6) with error term 0(A^~), as a consequence of an integral formula for M(A). For 
the Heisenberg subLaplacian we obtain 

M(A) = c'„A"+^ (1.8) 

for a constant given by an infinite series {c'l = ^5 C2 = 5^)- Here the Heisenberg 
group has Euclidean dimension 2n + l but Hausdorff dimension 2n + 2, so n + 1 is again 
half the Hausdorff dimension. For the Laplacian on a homogeneous tree of degree g + 1, 
we obtain an integral formula for M(A) and an asymptotic formula 

M(Ao + e) = c;'e=^/2 _^ ^(^3/2^ as e ^ 0+ (1.9) 

where \q = q + 1 — 2y/g is the bottom of the spectrum (note that the exponent | is 
independent of q). 

These results are all simple consequences of well-known formulas for harmonic analysis 
on these spaces, but our results place all these disparate facts together in a coherent 
context. 

In Section 4 we discuss some inhomogeneous examples. The simplest of these are the 
half-line and the Euclidean or hyperbolic half-space, with either Dirichlet or Neumann 
boundary conditions. In these examples K\{x,x) can be written as a sum of the kernel 
for the full space plus a term that decays as x — )■ 00. The decaying term thus contributes 
zero toward the average value, so M(A) is the same for the half-space as for the full 
space. More interesting examples are provided by unbounded open domains in R"". We 
are not able to obtain the desired results for such examples, but we sketch three different 
approaches to studying the problem. The issues here are quite technical since it is 
necessary to impose some hypothesis on the domain. This is an area that is certainly 
ripe for future development. 

In Section 5 we come to our second main heuristic idea, which we call the asymptotic 
splitting law. Suppose the space X and the Laplacian have a finite symmetry group G. 



4 



Then G has a finite number of inequivalent unitary representations {tTj } with dimension 
{dj}, and the dimension formula says 

^4 = #G. (1.10) 

i 

We may sort the eigenspaces into subspaces that transform according to the representa- 
tions {vTj}, and write M(A) as the sum of Mj{\) where each Mj{\) gives the contribution 
to the spectral mass arising from the representation tTj. The asymptotic splitting law 
heuristic is 

MAX) cR 

A"i M(A) #G ^ ' 

This statement is not always valid, as we demonstrate by simple counter-examples . For 
it to have any hope of being valid, we must have that most orbits under the G action on 
X be full orbits. But even more than this, we need the existence of a fundamental domain 
for the action that has boundary that is relatively small. We give a precise estimate for 
Laplacians on finite graphs for the difference between Mj[\)/M[\) and d'j/^G for all 
values of A. In some noncompact settings we can show that these quantities are exactly 
the same for all A. One interesting example we discuss is the Sierpinski gasket and related 
fractals. Here we are able to establish (1.11) with a much smaller error term than we 
would expect for nonfractal examples. 

In Section 6 we discuss similar questions when the symmetry group G is a compact Lie 
group. In this case there is no dimension formula, but we can still sort the eigenspaces 
according to the irreducible representations, and we can compare Mj{X)/Mk{X) with 
d^dl for different representations. 

In section 7 we discuss a noncompact fractal example, the generic infinite blowup 
of the Sierpinski gasket. The structure of the spectrum of the Laplacian here has been 
described by Teplyaev [Te] ; there is an orthonormal basis of compactly supported eigen- 
f unctions, but each eigenspace has infinite multiplicity, and the closure of the set of 
eigenvalues is a Cantor set. We present evidence that the spectral mass function ex- 
ists and has the same asymptotics as the eigenvalue counting function on the Sierpinski 
gasket. 

Finally, in Section 8, we discuss the relationship between the spectral mass function 
on a space X and a covering space X oi X. In typical examples X is compact and X 
noncompact, and we want to pass from information about the spectral mass function 
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M(A) on X to information about M(A) (and hence N{\)) on X from the identity 



where x projects to x and V is the group of deck transformations. The term corre- 
sponding to 7 = e, the identity, contributes to M(A) exactly the value M(A). We could 
conclude that M(A) has the same asymptotics as M(A) if we could show that the sum 
of all the other terms on the right side of (1.12) has slower growth rate. We give some 
simple examples where this is the case, but with the caveat that the infinite series is 
only conditionally convergent. We discuss the example of hyperbolic manifolds of finite 
volume, where such conclusions would be very interesting, but the technical challenge 
of estimating (1.12) seems difficult. For the example of the compact quotients of the 
Heiscnberg group, we note that the desired estimates have been obtained by Taylor [T] 
using heat kernel methods. Such an approach might well be more effective for dealing 
with other examples discussed in this section. For a direct approach to these quotients, 
see FoUand [Fo]. 

There arc many earlier words that arc closely related to our ideas. Hormander [H] 
gives estimates for the asymptotics of K\{x^ x) in the case of an elliptic pseudodifFerential 
operator on a manifold, but the estimates are not uniform in x at infinity or near the 
boundary (but also see [I]), and so do not immediately yield information about M(A). 
That paper also gives references to earlier work on the topic. The paper of Cheeger, Gro- 
mov and Taylor [CGT] gives estimates for the kernels of many functions of the Laplacian 
based on wave equations methods. Wave equation methods were introduced by Levitan 
[L] in 1952 (see also [SV1],[SV2]). Many of the expressions for K\ in the examples in 
section 3 may be found in our earlier paper [SI] and in the book of Taylor [T]. 

This paper is woefully incomplete, and is offered in the spirit of adventure. We hope 
that adventurous readers will be inspired to tackle some of the conjectures and questions 
that are scattered throughout the text. 

2 Spectral mass 

Let Ex denote the spectral resolution of a nonnegative self-adjoint operator A on LF'{X, jj), 
so Ex is the projection operator onto the [0, A] portion of the spectrum of ^4. In many 




(1.12) 



7er 
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cases Ex is an integral operator with a continuous kernel Kx{x,y): 

Exfix)^ f Kx{x,y)f{y)dii{y). (2.1) 
Jx 

For example, if A has a discrete spectrum of eigenfunctions ipj with eigenvalues 

A^j = (2.2) 
{</7j} an orthonormal basis of L^, then if Ex has finite dimensional range we have 

Kx{x,y) = ^ ipj{x)ipj{y) (finite sum). (2.3) 

A,<A 

In this case, 

J Kx{x, x)dii{x) = N{y) = #{Aj < A} = dim Range Ex- (2.4) 

More generally, we would like to define a spectral mass function M(A) that measures 
the "size" of projection Ex- We propose taking 

M(A) = Average(i\:A(a;,a;)). (2.5) 

Of course this is not a precise definition, since we have to say what wc mean by 
"Average", and this will have to depend on the context. Typically the average will be 
defined by 

1 f 

lim / Kx{x,x)diJ,{x) (2.6) 



for some "reasonable" sequence of sets Bk increasing to X, where the hmit is independent 
of the choice of the reasonable sequence. We may have to settle for liminf and lim sup 
in place of lim in (2.6). Two situations where the meaning of Average is obvious: 

(i) if ^{x) < oo then the average is just 

— — ■ I Kx{x,x)dii{x), in which case (2.7) 
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(ii) if Kx{x, x) is constant in X, then the Average is this constant. 

If X is a metric space we would hope to choose to be the baUs of radius k centered 
at some point z, with (2.6) independent of the choice of z. 

The following known result provides some basic properties of Kx{x,x). We provide 
the simple proofs for the benefit of the reader. 

Lemma 2.1. a) Kx{x,x) is nonnegative (possibly +00). 

b) Kx{x,x) is monotone increasing in X. 

Proof, a) Since Ex is a projection 

Kx{x,z)= f Kx{x,y)Kx{y,z)di^{y) (2.9) 
Jx 

If we set X — z and observe that Kx{y,x) = Kx{x,y) because Kx is self-adjoint, (2.9) 
yields 

Kxix,x)^ [ \Kx{x,y)\^dM- (2-10) 
Jx 

b) Suppose A' < A. Because Ex is a spectral resolution, Ex — Ey and Ex' are 
orthogonal projections {Ex — Ey is the spectral projection onto the [A', A] portion of the 
spectrum) with kernels Kx — Ky and Ky- The orthogonality implies 

/ {Kx{x,y) - Ky{x,y))Ky{y,z)dfi{y) ^ 0, (2.11) 
Jx 

so the same reasoning that leads to (2.10) shows 

Kxix,x)-Ky{x,x)= [ \Kxix,x)-Ky{x,y\^dfi{y) (2.12) 

Jx 

□ 

The lemma shows that if we can make sense of (2.5), then the spectral mass will be 
nonnegative and monotone increasing in A. 

It is interesting to compare the spectral mass with a more familiar object, the heat 
kernel on the diagonal. Let ht{x, y) denote the kernel for the heat semigroup e*^. There is 



8 



a vast literature on estimates for heat kernel, both on and off the diagonal ([G]). Typical 
results are upper 

ht{x,x) <ci^t) (2.13) 

and lower 

ht{x,x)>c,^t) (2.14) 

on diagonal bounds, for specific (or sometimes generic) functions ^{t), often a negative 
power $(i) = for some 6 > 0. Define 

M+(A) = limsup p / Kx{x,x)dii{x), (2.15) 

If 

M (A) = lim inf . , / Kx{x,x)diJ,{x) 



and 



//+(i) = limsup—^ / ht{x,x)diJi{x), (2.16) 
H~ (t) — lim inf / ht{x,x)dii{x). 

k-^oo li[Bk) Jb, 



Since we know 



poo 

ht{x,x)= j te-^^Kx{x,x)d\ (2.17) 
Jo 

wc may draw various conclusions relating these quantities. For example, the upper 
estimate (2.13) implies 

H'^it) < ci$(t), (2.18) 

and similarly (2.14) implies 

H-{t) > C2^t). (2.19) 

On the other hand (2.17) implies 

Kx{x,x) < ehi/x{x,x) hence (2.20) 
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A X / Kx{x,x)dfi{x) < ^ f hi/x{x,x)dfi{x) (2.21) 

/^l-Dfcj JBk fJ'[^k) JBk 

and finally 

M±(A) < eH^{l/\) (2.22) 
So the upper estimate (2.13) implies 

M+(A) < eci$(l/A). (2.23) 

Of course, if we know power law limits for the heat kernel, we can use (2.17) and the 
Karamata Tauberian Theorem to obtain power limits for Kx{x,x). 

One of the main tools is estimating the eigenvalue counting function A'"(A) is the 
Rayleigh quotient 

n(u) = f^, (2.24) 

where £ is the energy associated to the Laplacian and the measure via 

- J{Au)vdfx (2.25) 

The minmax characterization of eigenvalues leads to 

N{X) = max{dimL : 7^|L < A}. (2.26) 
If A and A' are two Laplacians and 

n(u) < n'(u) for all u, (2.27) 
where TZ' is the Rayleigh quotient for A', then (2.16) implies 

A^(A) > A^'(A). (2.28) 

Typically we start with one Laplacian and construct two simpler Laplacians Ai and 
A2 for which 

7^l(^^) < n{u) < n2{u). (2.29) 

If we can show that Ai and A2 have the same spectral asymptotics, then we may conclude 
that A also has the same spectral asymptotics, 

In the general context, the Rayleigh quotient is well-defined by (2.24) and (2.25). 
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Conjecture 2.2. Suppose A and A' are two Laplacians for which (2.27) holds. Then 

M^(A) > (M^)'(A). (2.30) 

This conjecture would allow us to use Reyleigh quotient estimates to obtain spectral 
mass estimates. It would appear that it would be easy to prove this conjecture in con- 
texts where the integrated density of states definition fll.7p is valid. It is also an easy 
consequence of the next conjecture. 

Conjecture 2.3. Let L denote a closed subspace of L'^{fi), and KL{x,y) the kernel of 
the orthogonal projection onto L. Then 

M"*" = max I lim sup — — — / KL{x,x)dfi{x) : 71\l < \\ (2.31) 

M~ = max I liminf / Kl{x, x)dn{x) : 71\l < A 



with the maximum attained for L = range E, 



3 Homogeneous spaces 

In this section we discuss some examples in which there is a transitive group of symmetries 
actiong on X that preserves the operator A. In such cases the function Kx{x,x) is 
independent of x, so the Average in (2.5) is just this constant value. 

Example 3.1: The Euclidean Laplacian 

Let A = -A on M". Then 

Kx{x,y) = j^[ e^(^-^Hrfe (3.1) 

It is well-known how to express Kx{x,y) in terms of Bessel functions, but for our 
purposes it suffices to note that the integrand is identically one where x = y, so 

Ka(x,x) = c„A"/2 for (3.2) 
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where 6„ is the volume of the unit ball in M". The constant is the same one that appears 
in the Weyl asymptotic law, but in this case the asymptotic law is exact for all values of 
A. 

Example 3.2: Lattice Graph Laplacian 

Let X = Z", the n-dimensional cubic lattice, and let A be the negative of the graph 
Laplacian 

n 

- A/(m) = ^(2/(m) - /(m + e,) - /(m - a)) (3.4) 

i=l 

where denotes the standard basis elements in Z". (Note: it is sometimes conventional 
to divide by 2n in the definition, but this only changes the computations below by a 
constant.) The theory of Fourier series tells us that reasonable functions / may be 
written 

/(^) = / E (3.5) 

where C denotes the cube : < |} in R", and the functions m — )■ e^'^™'^ are 
eigenfunctions of —A with eigenvalue X]r=i ^(1 ~ cos27r,^j) = Yl'i=i 4sin^ n^i. This shows 
that the spectrum of —A is [0,4n] and the spectral projection operator Ex is obtained 
from (3.5) by restricting the integral to 

Cn |^4sin2 7rCi < a| . (3.6) 

Thus the spectral mass function is the volume of the set (3.6). When n = 1 this gives 
exactly 

M(A) = ^sin-(^^j. (3.7) 

Note that this is asymptotic to ^ as A — > 0. For higher values of n wc do not have an 
exact formula for M(A), but for small values of A we may approximate sin^ n^i by (tti^i)^, 
so 

M(A) ~ c„A"/2 as A ^ (3.8) 

for the same constant (3.3). In other words, the Weyl asymptotic law holds for the 
bottom of the spectrum. 
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Example 3.3: The Hyperbolic Laplacian 



Let A = —A on i/", the hyperbolic n-space, with metric d{x, y). It is known that the 
spectrum is {^^Y ' write 

Kx{x,y)^ I ' ^t{d{x,y))dt (3.9) 

for an exphcit spherical function (/^^ ([Ta]). For our purposes it suffices to know 

r((^)%.t) " 



9'i(0) 



Wn-l 

(27r)" 



r{it) 



(3.10) 



n— 3 
2 



71 — 4 

^ ttanhvrt t] 1 + \ j + - I 



(2 



n odd 



n even 



where w^-i denotes the surface measure of the unit sphere (so Wn-i = nBn). Thus the 
spectral mass is exactly 

M(A) = Cnu I 
Jo 

Asymptotically we have 



r{it) 



dt. 



(3.11) 



= CnA2 + 0(A 2 j as A ^ 00. 



(3.12) 



Note that the remainder term is smaller than expected. This may also be explained in 
terms of more refined Weyl asymptotic laws. For example, when n — 3 the exact formula 
is C3(A - 1)2. 

Example 3.4: The Heisenberg SubLaplacian 

Let Heis„ be C" x M endowed with the group law 



{z,t) o {z',t') = (^z + z',t + t' - ^Imz ■ z'^ 



(3.13) 
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with left invariant vector fields 



V - — - — 



(3.14) 



T 



d_ 
di 



and let A = for 

n 

£ = Y,iXj+y'), (3.15) 

the Heisenberg subLaplacian. The joint spectral resolution of the commuting operators 
£, iT is known, and we may write 



oo „A 
c — 1-1 I — n •'0 



(3.16) 



e=±l k=0 

where * denotes the group convolution. Here fx^k.e is a the spherical function 



<^A,ifc,e(^, t) 



+1 



exp 



—ieXt 



exp 



n + 2kj \A{n + 2k) 



-\\z\ 



T n—1 



X\z\ 



2{n + 2k) 



(3.17) 



and L? ^ are the Laguerre polynomials characterized by the generating function identity 



^r'=L^-^(s) = (1 - r^e-^ for |r| < 1. 

k=0 

Thus the spectral mass is 

oo „x 
c — ui I — n •'O 



(3.18) 



(3.19) 



e=±l fe=0 ■ 

oo „A 



= 2E 



(2t)"+"(» + 2*:) 



n+1 



Lri(o)ci.. 



Note that (3.18) implies 



1 / d 



r=0 



n + /c — 1 
A; 



(3.20) 
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and so 

It is clear that the infinite series converges since the terms are 0{k~'^). It appears hkely 
that one can evaluate the constant in (3.21) in terms of the values of the zeta function 
for even integers. For example, when n = 1 it is 



(27r) 

and for n = 2 it is 



^' + 1 _ 2 1 1 



3(27r)3 ^ (2 + 2A;)3 3(27r)3 8^'' 5767r 



Example 3.5: Laplacian on Homogeneous Trees 

Let Tg denote the homogeneous tree where each vertex has q + 1 distinct neighbors, 
with q > I, and let A = —A where 

A/(x) = 5](/(y)-/(x)). (3.24) 

yr^x 

A priori, the spectrum of A lies in [0, 2{q + 1)], but in fact it turns out to be the smaller 
interval [q + 1 — 2^, q + 1 + 2y/q]. The spectral resolution is due to Cartier [Ca] and 
may be found in [F-TN]. 

Define the c-function by 

c{z) 

and the spherical functions 

cj),{x) = c{z)q-'^^''^ + c(l - ^)g-(i-^W^) (3.26) 

where d{x) denotes the graph distance to some fixed reference vertex xq G T. The 
spherical function (3.26) is an eigenf unction of A with eigenvalue 

A = g + l-(g" + g^-"). (3.27) 



q + l J q 



,2-1 



for ^ e C 



(3.25) 
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It turns out that the only contribution to the spectral resolution comes from Rez — |, 
and we will write z — ^ + it ior < t < so (3.27) becomes 

A = g+ 1 - 2^cos(tlogg). (3.28) 



The full spectral resolution may be written 



log q 



(p^+it * f{x)dm{t) 



where the convolution is 



* f{x) = (l)i+i^{d{x, y))f{y) 



(3.29) 



(3.30) 



and 



"^"^^'^ ^ 2^^ 



c[ - + it 



dt = y + l)loggsin^(Uogg) 
- V. . . 27r((9-l)2 + 4gsin2(Uog?)) ^ ^ 

The spectral projection Ex is obtained by restricting the integral on the right side of 
(3.29) to the values of t with A < t, in other words 



< t < 



log? 



cos 



q +l- X 



(3.32) 



The definition of the c-function was chosen to have c{z) + c{z + 1) = 1 so that (^^(0) — 1, 
so 

Kx{x,x) = / dm{t) (3.33) 

where Ix is the interval (3.32), and this is the spectral mass function. Note that M(A) = 
for A < g + 1 — 2y/q (below the spectral gap), and just above the spectral gap A = 
q + 1 — 2y/g + e for small e. 



Jo 7r(g - ly 



3TT{q - 1)2 



^3/2 + 0(6^/2). 



(3.34) 
(3.35) 
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4 Inhomogeneous spaces 



In this section we discuss examples of spaces that are not homogeneous, forcing us to 
come to grips with the notion of "Average" in (2.5). 

Example 4.1: The half-line 

Let X be the half-line [0, oo] and ^4 be — ^ with cither Neumann or Dirichlet boundary 
conditions at the origin. Then Ex corresponds to Ex on the full line with the function 
extended by even or odd reflections. In other words, 

Kx{x, x) = Kx{x, x) ± Kx{x, -x) (4.1) 

where Kx is the kernel for the full line (Example 3.1 with n — 1). Explicitly, 

lsWA(x-y) 
TT X — y 



The obvious definition of "Average" here is 



1 ''^ 



M{\)= lim - / Kx{x,x)dx. (4.3) 
-R^oo R Jq 

Since Kx{x, —x) = ^ ^'"^^ tends to zero as x — >■ oo, it clearly does not contribute to 
(4.3) and so we obtain the same result ciVX as for the full line. 

Example 4.2: The half-space 

Let X = X [0, oo) and A be —A with either Neumann or Dirichlet conditions on 
the boundary. If we write x' — {xi, . . . , x„) e M" and < Xn+i < oo then the analog of 
(4.1) is 

Kx{{x',Xn+l),{x',Xn+l)) = Kx{{x' , X^+l) , {x' , Xn+l)) ± Kx{{x\ Xn+l), {x' , -Xn+l)) (4.4) 

where Kx denotes the kernel for R""*"^. It is well-known that Kx may be expressed in 
terms of Bessel functions and Kx{{x' ,Xn+i), {x', —Xn+i)) tends to zero as Xn+i — >■ oo. So 
if we define "Average" in the expected way, 

M(A) = lim — — ^— — / Kx{{x',Xn+i),{x',Xn+i))dx'dxn+i (4.5) 
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then Kx{{x',Xn+i), {x\ —Xn+ij) in (4.4) does not contribute to the hmit in (4.5), so we 
get the same answer c„+iA^~ as in the case of M"^^. There are many other ways we 
could define the "Average" that would lead to the same answer. On the other hand one 
could imagine some very perverse definitions that would take averages over regions that 
eventually exhaust the half-space but unduly weight points close to the boundary, and 
lead to different answers (or even nonexistence of limits) . 

Example 4.3: Disjoint Cubes 

Let {Cj} be an infinite sequence of open cubes in IR" whose closures are disjoint, 
let X — {}■ Cj and let A be the restriction of —A to each cube with either Neumann 
or Dirichlet boundary conditions on the boundary of each cube. This is not usually 
considered an interesting example since X is not connected, but it will enable us to 
develop some insight that may be carried over to more interesting examples. Let Sj 
denote the side length of Cj, so vol(Cj) = s^. 

To begin we make the following assumptions 

(i) Ylij -Sj = so X has infinite volume, but 

(n) there exists e > such that Ylfs <e^l < oo, so that all sufficiently small cubes 
make only a finite contribution to the volume. 

We define "Average" as the limit of averages over UjLi Cj as N ^ oo, but in fact 
this limit may not exist, so we define upper and lower averages 

1 ^ r 

M+ = limsup — ^ / Kx{x, x)dx (4.6) 

and M~(A) with limsup replaced by liminf. Now the projection operators E\ are given 
on each Cj as just the standard Neumann or Dirichlet Laplacian projection operators on 
Cj, so in fact 

/ Kx{x,x)dx = Nj{\) (4.7) 
where Nj{\) is the eigenvalue counting function on Cj. 
We know that 

AT, (A) = c„s^"A"/2 + ej{\)s';r'^ (4.8) 

where the error term ej{\) tends to zero as A — )■ oo. However, the error estimate is not 
uniform across all cubes. The two observations we need are that in any cases the error 
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term ej{X) is uniformly bounded for all cubes, and the vanishing ej{X) as A — > oo is 
uniform on all cubes with Sj > e (for example, for any fixed A, if Sj is small enough then 
Nj[X) = 1 (Neumann) or (Dirichlet)). The first observation together with assumptions 
(i) and (ii) mean that the small cubes with Sj < e do not contribute to the limit defining 
M^(A). The second observation then yields estimates 

c„A"/' - e-(A)A"/' < M-(A) < M+(A) < c„A"/' + e+(A)A"/' (4.9) 

where — )■ as A ^ oo. Thus (4.9) is our version of Weyl asymptotics in this case. 
(Better estimates are possible, but we will not discuss them here.) 

On the other hand, if we violate assumption (ii) then we can end up with essentially 
meaningless results. Suppose, for example, that (i) holds but limj_^oo Sj = 0. Then for a 
fixed A we will have all but a finite number of cubes too small to support nonconstant 
eigenfunctions with eigenvalues < A. Thus Nj{X) = for Dirichlet boundary conditions 
for all but a finite number of j's, so 

M+(A) < hmsup — ^ = 0. (4.10) 

For Neumann boundary conditions with Nj{X) = 1 for all but a finite number of j's we 
have 



M-(A) > hminf-^^ = +oo. (4.11) 



Example 4.4: Open subsets of R"' 

Let X be an unbounded open subset of R" with smooth boundary, and consider 
—A with either Dirichlet or Neumann boundary conditions. Under some additional 
assumptions we expect to obtain the estimate (4.9). This is weaker than the result for 

the half-space in Example 4.2, so again we could hope for the stronger conclusion under 
more stringent hypotheses. This is likely to be a long term project, so we restrict ourself 
here to outlining three prospective approaches to the problem: 

(a) If Conjecture 2.3 is valid, then we can try to find unions of cubes bracketing X, 
namely 

IJQCXCIJC;. (4.12) 

j j 

and use the estimate (4.9) for Example 4.3. For the inner approximation we can arrange 
for the closures of Cj to be disjoint, but not for outer approximation, so this would 
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require a minor adjustment of the argument (but this is no different from the case of 
bounded domains). The main drawback to this approach is verifying condition (ii). For 
example, if X is the region above the graph of a smooth function, then we would need 
to assume that the function rapidly approaches a constant at infinity to guarantee that 
small cubes have finite total volume. It is likely that a more careful argument could 
allow a weakening of condition (ii) to allow a more reasonable class of domains. Since 
Conjecture 2.3 is needed to pass from elementary Rayleigh quotient estimates that follow 
from (4.12) to estimate for spectral mass, it is not worth working out the details of this 
approach until the status of the conjecture is resolved. 

(b) Heat kernel estimates of the form (2.13) and (2.14) are known to hold, so we can 
use these to obtain estimates (4.9). The only technical problem is to find the hypothesis 
on X that will make the constants in (2.13) and (2.14) uniform across the unbounded 
set X. This approach is too cride to yield the more precise asymptotics that we found 
for the halfspace. 

(c) Wave equation methods, as described in [CGT], [L], [SVl], [SV2] can be used 
to compare functions of the Laplacian on the ambient Euclidean space. There are two 
technical issues here. The first is that we cannot use the spectral projection operator 
directly because the sharp cutoff at frequency A is implemented by ■0(— A), where is a 
discontinuous function whose Fourier transform decays too slowly at infinity. So we have 
to bracket E\ by 

V'i(-A)<£;a<V'2(-A) (4.13) 
where ■^i and ■02 are soft cutoff functions in T>. We may take to be even functions so 

i)j{s)^ / 7p{t) cos{st)dt and (4.14) 
Jo 

POO 

^j{-A)^ ifj{t)cos{t^/^)dt (4.15) 
Jo 

represents ipj{—A) as an integral of wave propagator operators cos{t\/—A). Because of 
the finite propagation speed of cos{t^/—A) we have the equality of the kernel cos(t\/— A) 
for the Laplacians A on X and A on as long as we stay a distance of t away from 
the boundary. Let us denote by K^^ the kernel of the operator ■0j(— A), and by K^^ 
the kernel of the operator ■0j(— A), and by Wt and Wt the kernels of cos(i-\/— A) and 
cos{tV^). Then (4.13) gives 

K^x\x,x) < Kx{x,x) < K^^\x,x) (4.16) 
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while (4.15) gives 



(4.17) 



y=x 



K^\x,x)= / ^^{t)Wt{x,y)dt 



(4.18) 



y=x 



(since the kernels Wt{x^y) are singular at y = x we must first integrate before setting 
y = x). The finite propagation speed implies 



where d{x) denotes the distance to the boundary. Since ipj vanishes rapidly at infinity, it 
should be possible to find a favorable estimate for the right side of (4.19), under suitable 
assumptions on X, using the method of [CGT]. This remains a technical issue. Since 
k'^^\x^x) is independent of x and has the desired 0(A"/^) asymptotics, (4.16) would 
then yield (4.9). Once again there is no hope of getting the existence of M(A) form the 
estimates (4.16). 

Example 4.15: The hyperbolic half-space 

Consider half of iJ" (Example 3.3) with either Dirichlet or Neumann boundary condi- 
tions (in the Poincare ball model of this would be a half ball). By the same reasoning 
as in Example 4.2 we will have exactly the same asymptotics, in this case (3.12) as for 
the full hyperbolic space. 

More generally we could consider other unbounded domains in using methods 
(b) and (c) sketched in Example 4.4. Even more generally we could look at Riemannian 
manifolds with positive injectivity radius and positive curvature bounds. It is not clear 
if we would require two-sided or only one-sided curvature bounds, and whether Ricci 
curvature bounds would suffice. 

5 Finite symmetry groups 

Suppose the space X has a finite group G of symmetries (x — )■ gx for each g E G) that 
preserves the measure /j, and the operator A {Af{gx) = A{fg){x) for fg{x) = f{gx)). 




y=x 



(4.19) 
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Let {ttj} be a complete set of inequivalent irreducible representations of dimensions dj. 
We know the dimension formula says Ylj^^'j — H^G. Wc may also split the span LP'{X) 
into an orthogonal direct sum ©jL|(X) where consists of functions that transform 
according to the representation tTj. In other words, / G if {/gjgGG spans a space on 
which the action of G is equivalent to tt^. We may then define spectral mass functions 
Mj(A) by restricting the spectral projections E\ to the spaces Lj{X). Clearly 

M(A) = ^M,(A), (5.1) 

i 

and a natural question that arises is what can be said about the ratios Mj(A)/M(A). A 
simple answer suggests itself: 

Asymptotic splitting law heuristic 

This heuristic is certainly not universally true. An obvious necessary condition is 
that most orbits under the action of G be full orbits (have the same cardinality as G) 
since orbits that arc not full will not support functions that transform according to all 
irreducible representations of G (in other words, some of the spaces L| will be too small). 
But even this condition is not sufficient, as the following simple example shows. Suppose 
X is the 3-regular graph consisting of two copies of Z^v, {xo,Xi, . . . ,xn_i,xn = Xq} 
and {yo,yi,---,yiv-i,yjv = Vo} with xj ~ Xj-i, Xj ~ Xj+i, yj ~ %_i, yj ~ yj+i and 
Xj ~ yj. Let G be the two element group consisting of the identity and the permutation 
Xj < — > yj for all j. Clearly the two spaces Ll{X) and Lf^X) consist of the even functions 
f{xj) — f{yj) and the odd functions f{xj) — —f{yo)- If ^ is the standard negative 
Laplacian 

(Af{xj) = 3f{xj) - f{xj_i) - f{xj+i) - f{yj) 
[Afiy,) = 3/(1/,) - f{y,.,) - f{y,^,) - f{x,) 



then it is easy to see that 



f(xA = eS-^'^/^ ^ ^ 

, 5.4 



are the even (+) and odd (-) eigenfunctions with eigenvalues 




(5.5) 
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(the multiplicites are 2 unless k — or k — y)- Note that the even spectrum hes in 
[0, 4] and the odd spectru, hes in [2, 6]. In other words, the bottom third of the spectrum 
consists entirely of even spectrum, the top third consists entirely of odd spectrum, and 
only in the middle have do the two interweave. Although this is a finite example, so the 
limit in (5.12) is meaningless, it clearly violates the spirit of the heuristic and may be 
used to construct actual counterexamples. 

It is easy to sec why this example misbehaves, and then to formulate conditions that 
rule out this sort of problem. If we split X into two natural fundamental domains for the 
action of G. namely the two copies of Zat, there are too many edges connecting them. 
What we want instead is a subdivision into fundamental domains with relatively few 
edges connecting them. We consider first the case of Laplacians on finite graphs. 

We say that F C X is a fundamental domain for the G-action on X if 

X=\JgF (5.6) 

g€G 

We do not insist that the images gF be disjoint, but of course we want the overlaps to 
be small. We define the boundary dF to consist of all points in F that either lie in, or 
have an edge connecting to some gF with g not the identity. 

Note that #(9P may vary depending on the choice of the fundamental domain, we 
would like to make a choice that gives close to the smallest possible value. 

Lemma 5.1. Let A be any G-invariant Laplacian on L'^{X,ii) for a finite graph that is 
nonnegative and self-adjoint, specifically 



for some positive conductances Cxy. Then 



ci? d'^.ifdF 



Proof. We may work with A'"(A) and Nj(X) (defined in the obvious way) since M(A) = 
N{X)/ ii{X), etc. Form the minmax characterization of eigenvalues we have 

iV(A) = max{dimL : Teii < A} (5.9) 
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where 



n{u) = ^^^^^ (5.10) 

2_^\u{x)\^^{x) 

X 

denotes the Rayleigh quotient and L denotes a subspace of LF'iX). Similarly 

Nx{y) = max{dimL : L C L]{X) and 7^|L < A}. (5.11) 
Let Vq denote the Dirichlet domain on F, namely 

Valine L'^{F) : u\dF = 0} (5.12) 

and let 

N^^\X) = max{dimL : L C Pq and 7^U < A}. (5.13) 

The key observation is that for each u E T>o and each representation ttj, there are 
linearly independent extensions of u to X and the Rayleigh quotient of each extension is 
the same as Tl{u). This immediately yields the estimates 

4Ar(o)(A) < Nj{X) < 4(ArW(A) + ^dF) (5.14) 

and summing over j we obtain 

#G'iV(°)(A) < N{X) < #G'(Ar(o)(A) + ifdF). (5.15) 

Combining (5.14) and (5.15) yields 



Nj(X) 



< dMF, (5.16) 



and (5.8) follows from (5.16) by dividing by Ijl{X). □ 

Although the lemma only deals with finite graphs, it does show that it is not necessary 
to go to very large values of A in order to get desired splitting to hold with small error; 
all that is needed is that be small in comparison to M(A). 

Example 5.1: Sierpinski gasket 

The Sierpinski gasket {SG) fractal in an equilateral triangle has the dihedral symme- 
try group ^3 of the triangle acting on it, and the group preserves the Kigami Laplacian 
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([Ki],[S2]). Since the spectrum of this Laplacian is known in great detail [GRS] it is 
possible to verify the asymptotic splitting law heuristic directly. However, Lemma 5.1 
together with the method of spectral decimation ([FS]) gives the result with a very small 
error estimate. SG is the invariant set for the iterated function system (IFS) consisting 
if the three contraction mappings Fi{x) = + ^Qi, where {qi} are the vertices of the 
equilateral triangle, 

SG^[jF,{SG). (5.17) 

i 

We construct a sequence of graphs that approximate SG as follows: Fq is the complete 
graph on the vertices {Qi}, and inductively 

F^ = Uf,(F^_i). (5.18) 

i 

Treating {qi} as the boundary of SG and each graph, wc note that all the other verices 
of F^ have degree 4, and so we may define a graph Laplacian 

Amu(x) = J2^u{y) - u(x)) (5.19) 

m 

at nonboundary points x with either Dirichlet or Neumann boundary conditions, and this 
is self-adjoint for the counting measure /i. The Kigami Laplacian on SG is a renormalized 
limit of these group Laplacians, 

3 

Au^- lim ^"'AmU (5.20) 

2 m— >-oo 



The method of spectral decimation says that the spectrum of A is a limit of the spectra 
of Ajn in a very precise sense: for any bottom segment [0, A] of the spectrum of A there 
exists m ^ j such that the eigenfunctions restricted to F^ are eigenfunctions of A^ 
and fill out the bottom segment [0, A'] of the spectrum of A^ (here A' ^ S^A). Since the 
action of the dihedral group preserves all these structures, wc have A^(A) = N^"^^{X') 
and Nj{X) = Nj"^\X') for the three representations ttq = trivial representation with do = 
1, TTi = alternating representation with di = 1 and 112 = 2-dimensional representation 
with d2 = 2 (note 1^ + 12 + 2^ = 6 = i^D^). So the ratio Nj{X)/N{\) on SG is exactly 
equal to the ratio N^^\X') / N^"^\X) on the graph T^n, and we may use (5.16) to estimate 
the difference from dj/6. We note that the natural fundamental domain F (either for 
SG or the graphs) is just 1/6 of the triangle bounded by the perpendicular bisectors of 
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the triangle. Now it is easy to see that #9F = m for which is 0(log A). Thus (5.16) 
imphes 

7V.(A) = iiV(A) + 0(logA). (5.21) 
In this case N{X) = 0(Ai°g3/iog5)^_ 

Another way to explain the small reminder in this case is that dF on SG is just a 
countable set, hence a set of dimension zero. 

In fact it is possible to improve the error estimate in (5.21) for the 2-dimensional 
representation 712 so that the error is zero! Out of the first 3A^ eigenf unctions, exactly 2A'" 
are associated to 772 (this requires a proper sorting of the high multiplicity eigenspaces) . 
Another way of saying this is that the eigenfunctions come in triples, two associated with 
772 and the remaining are associated to either ttq or tti. This can be seen by examining 
the explicit description of the spectrum via the method of spectral decimation [GRS], or 
by a Rayleigh quotient argument given in [ASST] for a related fractal, the pentagasket 
(spectral decimation does not apply to this fractal, but it is post-critically finite (PCF)). 
In fact, the reasoning in [ASST] yields the following. 

Proposition 5.2. For a large class of PCF fractals with dihedral Dk symmetry group, it 
is possible to sort the eigenfunctions (with increasing eigenvalue order) into groups of k, 
with each group containing all [^^] 2-dimensional representations and half (one when k 
is odd, and 2 when k is even) the 1-dimensional representations. 

In particular, for the 2-dimensional representations 

A^,-(A) = ^Ar(A) ± k, 
while for the 1-dimensional representations we only have 

NA>) = ^N{\)-rO{\og\). 

Even more surprising, there is experimental evidence that for the Sierpinski car- 
pet with D^^ symmetry, something very close to Proposition 5.2 holds: there are some 
"mistakes" in the groupings but the mistakes eventually get "corrected". This will be 
discussed in [BKS]. 

Example 5.2: The 2-torus 
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The dihedral group D4 acts on the 2-torus / 1} and is a symmetry group of the 
standard Laplacian. It has four 1-dimcnsional representations and one 2-dimensional 
representation (1^ + 1^ + 1^ + 1^ + 2^ = 8). The eigenfunctions of the Laplacian are 
just e^'r*('^i^i+^2a;2) -with (ki,k2) E Z^, and the action of D4 on 7? describes the action 
on the eigenspaces. Note that a generic point in 7? has an eight element orbit, and 
the associated 8-dimensional eigenspace (it may be a subspace of a larger eigenspace) 
sphts exactly as expected. On the other hand, a four element orbit corresponds to 
just the trivial representation. In other words, if we consider the fundamental domain 
F — {{ki,k2) ki > k2 > 0} with boundary dF — {{ki,k2) : ki — k2 > or k2 — 0}, 
then 



< # k2) edF:kl + kl < = 0{\"^). (5.22) 



The same estimate is valid on the square with either Dirichlet or Neumann boundary 
values. 

It seems likely that this example is typical for Laplacians on compact manifolds. 

Conjecture 5.3. Consider the Laplace- Beltrami operator on a compact Riemannian 
manifold of dimension n with a finite symmetry group G (if the manifold has a boundary, 
assume the boundary is sufficiently regular, and take Dirichlet or Neumann boundary 
conditions). Then 

Nj{X) - -j^N{X) = 0(aV) 05 a ^ 00. (5.23) 
wG 

It should be possible to prove this conjecture either by modifying a proof of the Weyl 
asymptotics to take into account the G action, or else to realize the spaces L'j{X) as 
spaces on a fundamental domain F with certain boundary conditions, and then observing 
that the Weyl asymptotic law holds for those boundary conditions ([SV1],[SV2]). 

What is perhaps more interesting is that in the noncompact setting we may be able 
to strengthen the statement to their exact identity 

d^- 

Mj{X) = -j^M{X) for all A. (5.24) 

#G 

We illustrate this in the following examples. 
Example 5.3: The squEire lattice 



27 



We return to Example 3.2, taking n = 2, and consider the D4 symmetry. In this 
case the parameter space C for the cigcnfTinctions is the square [— |, |] x [~^) |] 

acts also on C. Almost every orbit is a full eight element orbit, and so all repre- 
sentations are involved in the 8-dimensional eigenspace spanned by e^'^*(='=™'i^i='=™26) Q^Yid 
g27r«(±m25i±mi52)^ Yqx example, the trivial representation ttq occurs in the average of these 
eight eigenfunctions, which is | cos 27rmi(^i cos 27rm2C2 + | cos27rm2Ci cos27rmi(^2. Thus 
Mo (A) is the average over 7? of 

j ^ (cos 27rmi^i cos 27rm2^2 + cos 27rm2^i cos 27rmi^2)^o?C- (5.25) 

But the average of the integrand is | for almost every and it is straightforward to inter- 
change the average and the integral to obtain Mo(A) = |Af(A). The other 1-dimensional 
representations are given by similar expressions, so (5.24) holds for all of them, and we 
obtain the same result for the 2-dimensional representation by subtraction. 

Example 5.4: The plane with dihedral symmetry 

We return to Example 3.1 with n — 2 and the dihedral symmetry group Djv. For the 
trivial representation ttq 

KT{^^y) = ^ E ^A(:r,72/) (5.26) 

where Kx is given by (3.1). The contribution to K^^\x^x) from 7 = identity is exactly 
-^K\[x,x). For all other 7, \x — ^x\ on average is large, and Kx{x,^x) tends to zero 
as \x — 7x1 — >■ 00 because of the decay of Bessel functions Jq. Thus Mq{\) — ^M{\) 
exactly, and a similar argument establishes (5.24) for all representations. 



6 Lie group symmetries 

Suppose the group of symmetries of G considered in the previous section is a compact Lie 
group. Then the family {ttj} of inequivalent irreducible representations of G is infinite, 
but the dimensions {dj} are all finite. So the definitions of L'j{x) and Mj{X) make sense, 
and (5.1) holds. But there is no dimension formula, and the ratio in (5.2) is likely to 
tend to zero, and so gives no information. Our heuristic in this case is 

^ -k as A ^ 00 (6.1) 



Mfc(A) d 
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with the hmit being uniform if we restrict j,k < N for any finite N. 

For this to have any hope of being vahd, we must have most orbits under G to 
be full orbits, or equivalently, the stabilizer subgroups of most points must be trivial. 
For example, the n-sphere S'^ has symmetry group SO{n + 1) (or 0{n + 1)), and the 
eigenfunctions of the invariant Laplacian are the spaces of spherical harmonics. When 
n > 3 there are representations of SO{n + 1) that do not even appear in L^(5'"), and in 
the case n = 2, even though all representations appear, the ratios in (6.1) tend to dj/d^. 
Of course S"^ — SO{n + 1)/ SO{n), so the stabilizer subgroups are nontrivial. 

li X — G and the operator is the bi-invariant Laplacian (Casimir operator), then 
(6.1) is trivially true, since the eigenfunctions are just the entry functions of tt^, and so 
have dimension exactly once A is larger than the eigenvalues associated to tTj and 

We can modify this example to get a more interesting class of examples by considering 
X — G X M where M is a compact Riemannian manifold and the operator is the sum 
of the Laplacians on each factor. Write {A^*^^} for the eigenvalues associated with the 
representations {ttj} and {A^^^} for the eigenvalues (repeated in case of multiplicity) of 
the Laplacian on M. Then 



Since N^^^ satisfies Weyl's asymptotic law, it is clear that N^^\X-a)/N^^\X-b) 1 
as A — >■ oo, so (6.2) implies (6.1). 

Similarly, we can consider X — G x M where M is a space with a Laplacian for which 
M^^^ (A) is well defined and satisfies a power law asymptotic behavior. Then 



and (6.1) holds. 

We expect similar results to hold for various warperd product metrics on G x M, 
provided the orbits under G remain uniformly bounded in size. The following example 
shows that we may need to modify the notion of Average in some cases. 

Example 6.1: Let X = with the standard action of the rotation group SO (2). 
The one-dimensional representations {ttj} associated to e*-'^ are naturally indexed by 



iV,(A) = ({#A;:A(^)<A-Af })d^^ 
= 7VW(A-AfV|- 



(6.2) 



(6.3) 
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j e Z, and L| consists of functions f{rcosv,rsmv) — e^^'"fo{r), and we project L'^{X) 
onto ^{X) by 



j=-oo 

The kernel of Ex then sphts into 



J = -0O 



(27r 

Using the identity 



(6.4) 



1 Z"^" 

— / f {r cos{v — u),r sm{v — u))e^-'^du — II j f{r COS v,r sin v) 
and we recover / from its projection by 

oo 

f{rcosv,rsmv)— 11^ /(rcos rsini>). (6.5) 



njii'A((?'cosv,rsini>),y) (6.6) 



and so we would hke Mj{X) to be some sort of Average of K)^'{x, x) given by 

1 

— / K x{{r cos{v — u),rsm{v — u))e^^^du. (6.7) 
27r Jo 

By (3.1) this is 

-1 pZTV 

/ e^"''"°^("-"-^)e-^"'''^°^("-^)e'^"sd^dsdii. (6.8) 

) Jo Jo Jo 



1 r 

— / e^*"°^("-")e'^"dii = iV,(i)e^^» (6.9) 
27r Jo 

twice, we transform (6.8) into 

lo^ Jo ^'■^^■(^^)^''^'~''^^~''''°'^''~^^^^^^* (6.10) 

and then into 

- J,{srfsds = ^ I Jjisfsds. (6.11) 

Since (6.10) is O (i) as r — >■ oo we cannot use the usual average on to define Mj{X). 
However, if define 

Mj(A)= lim 4 / K^^\x,x)dx (6.12) 
R-^oo R J\^\<ji 

then the limit exists and equals -\/A/27r^ for all j (see [S] section 3). Thus (6.1) holds. 
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7 A fractal example 



In this section we discuss an example where the operator is a fractal Laplacian. For 
compact fractals, where the Laplacian has discrete spectrum, it is usually the case that the 
eigenvalue counting function has a power order growth, A^(A) ~ A", but the ratio A^(A) / A" 
usually does not have a limit [KL]. For example, for the Kigami Laplacian ([Ki], [S2]) on 
the Sierpinski gasket (SG), with either Dirichlet or Neumann boundary conditions, the 
power a = log 3/ log 5, and the ratio approaches a discontinuous multiplicative periodic 
function as A ^ oo. The power is rather easy to explain, since a contraction of 1/3 
in measure increases the eigenvalue by a factor of 5. The value 2a has been mistakenly 
interpreted as a "spectral dimension" under the erroneous assumption that the Laplacian 
is an operator of order 2, whereas in fact it is an operator of order d + 1, where d — 
log3/log(5/3) is the dimension of SG in the effective resistance metric [S2]. 

The example we look at is a generic infinite blowup SG^o of SG [S]. It was shown 
by Teplyaev [Te] that the appropriately extended Laplacian Aoo on SG^o has pure point 
spectrum, and each eigenvalue has infinite multiplicity. More specifically, let A = {A^} be 
the discrete spectrum of the Neumann Laplacian on SG. Then 5"*A C A for all positive 

integers m. Let A^o = Um=o 5~™A. Then each A G A^o is an eigenvalue of — A^o with an 
infinite dimensional eigenspace Sx of functions (in fact there is a basis Sx consisting of 
compactly supported functions). The actual spectrum of Aqo is the closure of Aqo, and 
is topologically a Cantor set, but in fact the eigenspaces £x for A e Aoo give the spectral 
resolution (2.1) as 

Exf{x)^Yl / Kx'{x,y)f{y)dpi{y) for (7.1) 

Kx'ix,y) = ^(Pj{x)(pj{y) (7.2) 
j 

where {fj} is an ortonormal basis of Sy- 

For the Laplacian on SG with Neumann (or Dirichlet) conditions, the analogous sums 
(7.1) and (7.2) are finite, and it is known that 

M(A) ~ A"'0(A) as A ^ oo (7.3) 

where ip is a. discontinuous bounded function (also bounded away form zero) satisfying 
the multiplicative periodic condition iplBX) = ip{X). 
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To understand the spectral mass function for SGoo we need to simultaneously understand 
the averaging in (2.5) and the behavior of the infinite series in (7.1) 

The infinite blowup SG^o is the increasing union of finite blowups F~^o- ■ -oF"^ {SO) = 
SGw where w — {wi, . . . ,Wm) is a word of length \w\ = m with Wj G {0,1,2} and 
{Fo,Fi,F2} is the iterated function system (IFS) defining SG. The blowup is generic 
if the infinite word (^1,^2, . . .) is not eventually constant. The measure of SGyj is S'"'', 
and the average on SG^o is 

Ave(/) = lim 3" 

|«)|->0O 

if the limit exists. If we don't know that the limit exists we may consider upper Ave"^(/) 
and lower Ave~(/) averages by taking lim sup and liminf on the right side of (7.4). 

Now fix a value of w with \w\ = m and fix an eigenvalue A' G Aqo- The basis of 
£\i given in [T] is not orthonormal, but it can be broken into three parts. One part 
consists of functions that vanish on SGyj. Although this is an infinite set, it will give 
zero contribution to 

3-l"'l I Kx{x,x)dn{x) (7.5) 

S Gyj 

The second part consists of functions supported in SGm- This is nonempty only when 
A' G 5~"^'Am for some m', and then it is finite but may be large. Its contribution to (7.5) 
may be understood in terms of the 5™' A' eigenspace on SG. The third part consists of 
the functions whose support intersects both SG^ and its complement. We will show that 
there are at most six functions in this part for each A'. 

Lemma 7.1. For fixed w and X' there exits an orthonormal basis of the eigenspace S^my 
consisting of functions supported in either SGu, or its complement, together with at most 
six eigenf unctions {^j} whose support intersects both SG^ and its complement. 

Proof. Consider the restrictions to the boundary points of SGy, of the functions in the 
eigenspace and their normal derivatives at these points. Since there are a total of six 
values in all, this is at most a six dimensional space. Choose a basis {^j} for this 
space. Then any eigenfunction may be written as a linear combination of {^Pj} plus an 
eigenfunction that vanishes together with its normal derivatives at the boundary points 
of SGyj. Such an eigenfunction can be split into a sum of eigenfunctions supported in 
SG^ and its complement. We apply Gram-Schmidt to each of the parts separately to 
obtain the desired orthonormal basis. □ 



fdn 



SG 



(7.4) 
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Note that the only time that there will exist nonzero eigenfunctions supported in SG^ 
will be when 5"* A' e A. 

The contribution to (7.5) corresponding to A' is 

3- ( Ar,(A') + / \^A^)?dt^{x)\ (7.6) 

where N2{X') denotes the number of orthonormal basis in the second part, and {^j} arc 
the orthonormalizations of the functions in the third part. But A'^2(A') is essentially (plus 
or minus 3) the multiplicity of 5"* A' in A. Thus 

3""^ J]iV2(A')~3-™ ArA(5"*A') = 3-"^Ma(5"^A) (7.7) 

A'<A 5'"A'<5'"A 

where A^a is the eigenvalue counting function for —A on SG and Ma is the spectral 
mass function. In view of (7.3) and the identity (5"*)" = 3"*, it follows that (7.7) tends 
to A"'0(A) as m — >■ oo. 

We would like to estimate the contributions of the finite sum in (7.6) to (7.5). Now 

IsG I'fjl'^^l^i^) — ^® estimate 

f \^j\'di,<3"^m\l. (7.8) 

'J SGyj 

Using the standard Sobolev type estimate 

<c8{u,uY''' (7.9) 



\u\ 



oo 



if u has a zero, we obtain 

P,||oo<c(A')^/^. (7.10) 
Combining (7.8) and (7.10) yields the estimate 

3""^E / \^A^)?dii{x)<c{)^f'' (7.11) 

j=l J SGw 

for the finite sum in (7.6). We then need to sum (7.11) over all A' < A. 

The spectrum A of SG can be broken into segments of on the order of 2" distinct 
eigenvalues on the order of 5"^, n = 1, 2, . . .. So there are about 2" eigenvalues in A such 
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that A' ^5"" and A' < A means n — m! < log A/ log 5. This leads to the estimate 

A<A' m' „<m/+i2gA 

— logs 

m' 

1 I log 2 X — ^ / 

= cA2+isi5^2™, 



which is useless because the series diverges. 

Nevertheless, it is likely that the sum of the right sides of (7.11) is much smaller than 
the sum of the left sides. The idea is that while we cannot expect to improve on the 
global estimate (7.10), we only need to bound on SG^^ which is only a very small 
portion of its support for large m! . If the values of Ipj were uniformly distributed over its 
support, then the left side of (7.11) would be dominated by 3~™, and in place of (7.12) 
we would have 



\l<r\ m' ^ , . loe2 m' ^ 



(7.13) 



A'<A m' „<„/+l£g2 

— log 



Note only is this finite, but it gives a lower order of growth than the main terms A"'?/' (A). 
While it is overly optimistic to expect such uniform distribution for any one particular 
eigenvalue, one could reasonably expect it on average in the sum over A'. 

Conjecture 7.2. On SGoo, M(A) is well-defined and satisfies the same asymptotics (7.3) 
as on SG. 



Of course upper estimates for the heat kernel on SGoo of the form (2.13) with $(t) = 
are known (see [BP] and [S]), so (2.23) implies the upper bound 

M+(A) < ct". (7.14) 

This is weaker than the conjecture, but gives supporting evidence. 

Note that the integrated density of states definition (11. 7p easily yields the desired es- 
timates, so another approach to proving the conjecture would be to prove the equivalence 
of f ll.7p and the spectral mass in this case. 
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8 Covering Spaces 



Suppose X is a covering space of X with covering map tt and covering group F, and the 
operator on X is the hft of the operator on X. Can we relate the spectral mass functions 
M(A) and M(A)? 



If Kx{x, y) is the kernel of the spectral projection on X, then 

7er 



(8.1) 



is the kernel of the spectral projection on X. Whatever notion of average we have on X 
lifts to X. The summand corresponding to 7 = identity e in F will yield Af (A) when we 
average over tt~^{x) = 7r~^{y). So we may conclude that M(A) and M(A) have the same 
asymptotic behavior if we can show that 

Ave Kx('yz,z) 

has a slower growth rate than M(A). 



for 7 7^ e 



(8.2) 



Example 8.1: Let X be the circle of R / Z and X be circle R /A^ Z for some positive 
integer N, with 7r{x) — x mod 1. The eigenfunctions of on X are e^'^**^^ with /c e Z 
with eigenvalues (27rA;)^, so 



M(A) = A^(A) = 1 + 2 



Ta' 

27r 



(8.3) 



(where [x] denotes the greatest integer in x) 



On the other hand, the eigenfunctions of — 4^- on X are e^^"^^^/^ with A; G Z with 



eigenvalues {2nk/Ny, so 



M(A) = iiV(A) = l|l + 2 



nVx 

2n 



It is clear from inspection that 

M(A) = M(A) + 0(1) as A ^00 
but we would like to deduce this from (8.1) and (8.2). Here we have 



(8.5) 
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i'i<[¥; 



so (8.1) says 



(8.6) 



'^=° Ife|<[^] 

Of course, the 7 = term yields (8.4), while if 7 7^ the sum over k is bounded by N, 
so the contribution to (8.7) is uniformly bounded for all A. 

Example 8.2: Let X be again the circle M / Z but now let X = R. Then, as observed 
in Example 3.1, 

.i(x~v)-ej^ 1 sin ^/X{x - y) 



(27r) TT 



(x-y) 

and 



M(A) = — (8.9) 

TT 



SO (8.5) continues to hold. From (8.1) we have 



, , VX 1 sin y/X'y 
K^{x, x) = — + V (8.10) 

The infinite series in (8.10) is not absolutely convergent, but it is easy to see by summation 
by parts that it is uniformly bounded. 

Of course (8.10) is a consequence of the Poisson summation formula, and a similar 
description works for covering R"/ Z". The problem with nonabsolute convergence 
in (8.10) may be overcome by replacing the spectral projections with regularized approx- 
imations. 

Example 8.3: Compact hyperbolic manifolds 

Let X be the hyperbolic spaces iJ" of example 3.3, and let F be a cocompact discrete 
group of isometrics, so X is a compact manifold of constant negative curvature. Then 
Kx is given by (3.9), but now we need to have an expression for (pt{d{x, y)) for all values 
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of d{x,y), not just zero. The result is known and can be given exphcitly in terms of 
Legendre functions, but to simphfy the discussion we take n — 3, when 



and so 



Kx{x,y)^ / ipt{d{x,y))dt (8.12) 

1 sin \/A — 1 r — — 1 r cos \/X — 1 r 
27r^ sin /ir 

for ci(x, y) — cosh r. 
Note that in the sum 

X;,(x,x) = ^^A(7a:^,x) (8.13) 

the choice 7 = e corresponds to r = and gives the value C3(A — 1)^/^ for all x. Note 
that for 7 7^ e we have r 7^ 0, and (8.12) has growth of order 0(A^/^), and moreover the 
dependence on r has exponential decrease as r — )■ 00. The values of d{'-fx, x) will increase 
as 7 — )■ 00 in r, but r = cosh^^ d{'jx, x) will increase at a slower rate. The values also 
depend on x, but most likely this is not too significant. It seems plausible that 

J2Kx{^^,x)^ (8.14) 

El sin VA — 1 cosh~^ d{'-fx, x) — \/A — 1 cosh~^ (0^(70;, x)) cos a/A — 1 cosh~^(ci(7a;, x)) 

7^e 



27r2 (cosh ^ ci(7a;, a;))2y^(i(7a;, x^ - 1 



is of order 0(A^/^) independent of a;. Note that this may require further assumptions 
about the group F, and it may be that the infinite series is not absolutely convergent. 

Of course we know from the Weyl asymptotic law that the leading term of the asymp- 
totic behavior of M(A) is CsA^^^, with remainder term O(A^). The whole point of this 
approach is to reduce the remainder estimate to 0(A^/^). 

Example 8.4: Finite volume hyperbolic manifolds 

Consider the same setting as the previous example, but now do not assume that F 
is cocompact, but only that X — X /T has finite volume. We cannot assert that the 
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Laplacian has discrete spectrum, so we must use M(A) rather than N{X). However, since 
X has finite volume, the average in defining M(A) is simply 

M{X) = ^—I Kx{x,x)dfi{x). (8.15) 



Again, the 7 = e term in (8.13) contributes exactly C3(A — 1)^/^ to (8.15), so once again 
the problem is to estimate (8.14) by 0(A)^/^. Of course, even an estimate o(A^/^) would 
be interesting in this case, since we don't have a Weyl asymptotic law to fall back on. 

Example 8.5: Compact nilmanifolds 

Let X = Heis„ as in example 3.4, and let F be a cocompact discrete subgroup of 
Heis„. Note that we are dealing with the subLaplacian on X = X/T. For simplicity we 
will take F = Z^" x i Z in C" x M. Then 

^gZ2n e=-oo e=±l k=0 ^ ^ 

where (ps,k,€ is given by (3.17). 

The choice {(, £) = (0, 0) gives exactly M(A) given by (3.21), so it remains to estimate 
the contribution to (8.16) from all the other choices of (C,^)- From (3.17) we see this is 
explicitly the sum of 



^^^^0 (27r)-(n + 2fc)"+i 



X 2 cos f (--Im Cz] exp ( —-^^^^ L"" ^ ( ^ ds (8.17) 

\n + 2k\2 ^\-4{n + 2k)J ^ \2{n + 2k)JJ ^ ' 

over (C,^) 7^ (0, 0). It seems rather challenging to estimate this sum. 

In fact the desired estimates are obtained in [T] , p. 79-80 by first studying the heat 
kernel and then applying a Tauberian Theorem. The heat kernel has the advantage 
of rapid decay at infinity to avoid convergence problems. A direct approach to these 
quotients for certain choices of F is given in [Fo] . 

It would also be interesting to relate these results to the group representation decom- 
position of L^(Heis„/F) in [R]. 

Acknowledgements: The author is grateful to one of the referees who provided 
very detailed and helpful comments on the original version of this paper. 
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